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Abstract — All codes with minimum distance 8 and codimen- 
sion up to 14 and all codes with minimum distance 10 and 
codimension up to 18 are classified. Nonexistence of codes with 
parameters [33,18,8] and [33,14,10] is proved. This leads to 8 new 
exact bounds for binary linear codes. Primarily two algorithms 
considering the dual codes are used, namely extension of dual 
codes with a proper coordinate, and a fast algorithm for finding a 
maximum clique in a graph, which is modified to find a maximum 
set of vectors with the right dependency structure. 



I. Introduction 

Let F 2 denote the n-dimensional vector space over the field 
F 2 and let the inner product (•,•): F£ x FJ -> F 2 be defined 
in the natural way as (u,v) = Ya=i u i y i' where addition is 
in F2. The Hamming distance between two vectors of F 2 is 
defined as the number of coordinates in which they differ, 
and the weight wt(u) of a vector v € F 2 is the number of 
the nonzero coordinates of v. A linear binary [n, k, d) code 
C is a /s-dimensional subspace of F 2 with minimum distance 
d = min{wt(c) : c e C,c / 0}. A generator matrix G for 
an [n, k] = [n, k,d > 1] code is any matrix whose rows form 
a basis for the code. The orthogonal complement C of C 
in F 2 is called the dual code of C and is an [n, n — k, d- 1 ] 
code, where d 1 - is called the dual distance of C and n — k 
the codimension. A generator matrix H for the dual code is 
called a parity check matrix of the code C. In this paper we 
will say that a code C is an [n, k, d] d code if it is an [n, k, d] 
code with dual distance d , Further, two binary linear codes 
C\ and C 2 are said to be equivalent if there is a permutation 
of coordinates which sends C\ to C 2 . Throughout this paper 
all codes are assumed to be binary. 

A central problem in coding theory is that of optimizing one 
of the parameters n, k and d for given values of the other two. 
Usually this optimization is related to the following functions: 
n 2 (fc,<i) - the minimum length of linear codes for given 
minimum distance d and dimension k and d 2 (n, k) the largest 
value of d for which a binary [n, k,d] code exists. Codes 
with parameters [n 2 (fc, d), k, d] and [n, k, d 2 (n, k)] are called 
optimal. There are many reasons to study optimal codes. These 
codes are interesting not only for detection and correction of 
errors. Some of them have rich algebraic and combinatorial 
structure. The problems of optimality are strongly connected 
and can be considered as packing problem in statistics and in 
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finite projective spaces [ 10 1 . Unfortunately, all these problems 
are, as many others in coding theory computationally difficult 
0] and the exact values of n 2 (fc, d) for all d are known only 
for k < 8 [4 |. Tables with bounds and exact values for <i 2 (n, k) 
are given in [6| and 

Another application of optimal codes is directly related 
to the design method of cryptographic Boolean functions 
suggested by Kurosawa and Satoh 1111 . In this case the optimal 
linear codes have to be with largest possible dual distance. 
More precisely one have to study the function N(d,d ± ) as 
the minimal n such that there exists a linear binary code of 
length n with minimum distance d and dual distance d^. The 
investigation of N(d,d ± ) seems to be much harder than the 
investigation of n 2 (/c,d). There are some general bounds (see 
Jl2l) but these bounds can be reached only for a few values 
of d and d^. In a previous work, we studied N(d,d^) for 
d 1 - < d < 12 by computer using the package Q-EXTENSION 
0. With this package we attempted to construct generator 
matrices and classify codes with fixed parameters. Practically 
we had no success in the cases where the dual distance 
was more than 6. In the case of codes with fixed minimum 
distance larger than 2, it is quite natural to look at the duals 
of the codes with needed properties and to the parity check 
matrices. In other words, extending an [n — 1, k — 1, d] code 
C\ to an [n, k,d] code C 2 can be considered as extending an 
[n — 1, n — k] d to an [n, n — k] d code with one coordinate. This 
approach helps us to develop two different methods, which 
are much more effective when we study N(d,d ± ), and also 
n 2 (fc,d), for d = 8 and d = 10. For small dimensions, it is 
convenient to use a brute force algorithm that takes generator 
matrices of all inequivalent [n, k, d] d codes as input, extends 
them in all possible ways and then checks the constructed 
codes for minimum and dual distance and for equivalence. As 
a result we get all inequivalent [n + 1, k, > d] d codes. This 
method is of course impossible to use for larger dimensions k, 
because of the large number of possible extensions. To avoid 
this problem we use a second method for larger dimensions 
adopting a strategy for bounding of the search space similar to 
the strategy for finding a maximum clique in a graph suggested 
in 03). 

In this paper we present two algorithms which can be used 
for constructing of linear codes with fixed dual distance. We 
give classification results for all codes with minimum distance 
8 and codimension up to 14 and minimum distance 10 and 
codimension up to 18. We would like to refer to |7| and |8| 
for a detailed bibliography of works which study linear codes 
with minimum distance 8 and 10. 



2 



II. Preliminaries 

In this Section we give some properties of linear codes and 
the relations with their dual codes which help us to design the 
construction algorithms. 

The first proposition considers the even linear codes, namely 
the linear binary codes which consist only of even weight 
vectors. 

Proposition 1 If d > 2 is even and a linear [n, k, d] code 
exists, then there exists an even [n, k, d] code. 

Proof: If C is a linear [n, k, d] code, the code produced 
by puncturing C in one coordinate has parameters [n — 1, k, d 
or d— 1]. Adding a parity check bit to all codewords, we obtain 
an even [n, k, d] code. ■ 
Practically, we use the following corollary: 

Corollary 2 If d > 2 is even and even linear [n,k,d] codes 
do not exist, then no [n, k, d] code exists. 

Later on we give the definition and some properties of 
residual codes. 

Definition 1 The residual code Res(C,c) with respect to a 
codeword c € C is the restriction of C to the zero coordinates 
of c. 

A lower bound on the minimum distance of the residual 
code is given by 

Theorem 3 ([14 f Lemma 3.9) Suppose C is a binary [n, k, d] 
code and suppose c G C has weight w, where d > w/2. Then 
Res(C,c) is an [n— w, k — 1, d'] code with d' > d—w+\w/2~\, 

and on the dual distance by 

Proposition 4 Suppose C is a binary [n, k, d] code with dual 
distance d^, c G C, and the dimension of Res (C, c) is k — 1. 
Then the dual distance of Res(C,c) is at least d . 

There is also a well known elementary relationship between 
the minimum distance of a linear code and the parity check 
matrix. 

Proposition 5 A linear code has minimum distance d if and 
only if its parity check matrix has d linearly dependent 
columns but no set of d — 1 linearly dependent columns. 

The next proposition gives a connection between weights 
of the rows of generator matrices and columns of parity check 
matrices. 

Proposition 6 Any even linear code C has a parity check 
matrix whose columns have odd weights. 

Proof: Let G = [hP] be a generator matrix for C in 
standard form. Then every row of P has odd weight. And 
accordingly the parity check matrix H = [P T I n -k] has only 
odd weight columns. The sum of all rows of H gives the 
all-ones vector. ■ 



Proposition 7 If a linear code C has t codewords of weight 
1, then its parity check matrix has t zero-columns. 

Proof: Without loss of generality, let u = (1, 0, . . . , 0) £ 
C. Then Vv 6 C ± , = (u, v) = Y%=i u i v i — v i- ■ 
This means that if we have found all inequivalent codes with 
dual distance 8 and dimension 14, then we have also found all 
codes with minimum distance 8 and codimension up to 14. 
Another well known fact is that for d odd, 

| {inequivalent [n, k, d] codes} | 
> | {inequivalent even [n + 1, k, d+ 1] codes} | 

with equality only if the automorphism groups are transitive. 
We define for convenience the function L(k, d^). 

Definition 2 Let d > 3. Then L(k,d ± ) is the maximum 
length n such that a binary [n, k] d code exists. 

The next theorem holds. 

Theorem 8 Let C be a binary linear [n, k, d] d code. Then 
d>n-L{k- 1,^). 

Proof: Let C be a code with parameters [n,k,d] d and 
Cd be a codeword in C of weight d. If d < n — L(k — 1, d- 1 ) 
then the residual code Res(C,Cd) has dual distance at least 
d^, length n — d > L(k — 1, d^) and dimension k — 1. This 
is impossible because L(k,dj~) > L(fc,c?^) for d^ > d^- ■ 
In fact, the function L(k,d ± ) has already been investigated 
in another setting for dual distance greater than or equal to 4 
in connection with studies on K-caps in projective geometries 
since it is known that a K-cap in PG(k — l : q) is equivalent 
to a projective g-ary [n — k, k] d code with d > 4. And 
L(k,d ± ) have thus been considered in connection with the 
function 

/J,„{N, q) = the maximum value of k such that there exist 
a K-cap in PG(N,q), 

where q is the order of the underlying Galois field, in our 
case 2. More information on this can be found in the survey 
by Hirschfeld and Storme iflOll . 

III. Computational Algorithms 

We use mainly two algorithms for the extension of codes. 
If the dimension of the considered codes with given dual 
distance is small we can find the next column of the generator 
matrix of the new code relatively easy because we can 
represent any generator matrix G in the packing form as 
a vector Gb of n computer words. This algorithm, named 
BRUTEFORCE, can be described with the following steps: 

Algorithm BRUTEFORCE: 

INPUT: Ci n p - Set of all inequivalent [n, k] d codes 
represented by their generator matrices in packing form. 
OUTPUT: C out - Set of all inequivalent [n, k + l] d± codes, 
var a:array[1..2' c — 1] of integer; 
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In the beginning, C ou t is the empty set. For any code C r 
in Cm P with generator matrix in packing form G r do the 
following: 

1) Set a[i]:=l for any i. 

2) Find all linear combinations b of up to d — 2 column 
vectors of G r and set a[b]:=0. 

3) For all j such that a[j] = 1 extend G r with one coor- 
dinate, equal to j, to G' r . If there are no codes in C ou t 
equivalent to C' r (generated by G' r ) do C out := C ou tUC' r 

The big advantage of this algorithm is given by Step 2. In 
that step, all possible solutions for the (n + l)th column of 
the generator matrices are determined with approximately with 
53i=i (™) operations. Actually, to find all vector solutions for 
the (n + l)th column, we take all fc-dimensional vectors and 
delete those which are not solutions. We find all sums of less 
than d^ — 2 columns of the known part of the generator matrix. 
Each sum gives us one vector which is not a solution and have 
to be deleted. All remaining vectors are solutions. 

In Step 3, we use canonical representation of the objects. 
The main priority of the canonical representation is that the 
equivalence (isomorphism) test is reduced to check of coinci- 
dence of the canonical representations of the structures. In the 
case of many inequivalent codes, the computational time for 
comparing is growing fast. A technique for surmounting this 
problem is worked out. We split the set of inequivalent codes 
into a big amount of cells according to a proper invariant. 

To explain the next algorithm, we need the following 
definition. 

Definition 3 Let E be a set of k- dimensional vectors. 

1) We call E p-proper if all subsets of p vectors of E are 
linearly independent. 

2) Let M be a k x n matrix. The set E is called p-proper 
with respect to M if EL) {columns in M} is a p-proper 
set. 

Observe that, by Proposition [5] the columns of a parity check 
matrix for an [n,k,d] code form a (d — 1) -proper set. 

We consider the following problem: How to find a set of 
t = d — 1 binary vectors which have a certain property, i.e., 
(d — l)-proper subset of the set of all possible binary vectors 
with respect to a fixed generator matrix. To attempt to solve 
this problem in reasonable time, we adopt an idea suggested 
by Ostergard in [13] for finding a maximum clique in a graph 
in the algorithm EXTEND. 

Let C be an [n, k, d] d code with a generator matrix G in 
the form 



Let 



V* = {(l,x 2 ,x 3 , . . . ,x k ) : Xi 6 F 2 } 



G 



00 ... 



Res d {C) 



11. ..1 



X 



= [ A I X ] (1) 



where Resd{C) is a generator matrix of the residual [n—d, k — 
1, > code. Given that we know all such inequivalent 

generator matrices the problem is reduced to finding all (d 1 - — 
1) -proper sets X with respect to A of d— 1 binary vectors on 
the form (1, x%, ... , Xk) T ■ 



(Remark: If the dual code is even we may, by Proposition 
[6] reduce the search space to the set of odd-weight binary 
vectors). Delete from V* all linear combinations of d — 2, 
or less, vectors from A. The remaining set 



V = {vi,v 2 ,..-,v N } 



(2) 



is the search space for our search strategy. Now, for each 
integer 1 < i < N, let 



Vj = {vi,v i+ i, . . .,v N } 



(3) 



and let r be the A^-tuple, defined by r[i] = min{s, t}, 1 < i < 
N, where s is the size of the largest (d — l)-proper subset of 
Vi. First we consider (d — l)-proper subsets with respect to 
A of Vm that contain the vector vm, this obviously is {vn}, 
and we record the size of the largest proper subset found up 
to now in the tuple r, so r[iV] = 1. 

In the i-th step we consider (d — 1) -proper subsets with 
respect to 4 in V; containing Vi and record the minimum 
between size of the largest proper subset found up to now and 
t in r[i] (row *** in the algorithm). 

The tuple r for the already calculated steps enables the 
pruning strategy for the search. Since we are looking for a 
proper subset of size t = d — 1, and if the vector Vi is to 
be the (size)st vector in the subset and size + r[i] < t, then 
we can prune the search (row * in the algorithm). When the 
search terminates, the size s of the largest (d 1 - — l)-proper 
subset with respect to A of V or t (if t < s) will then be 
recorded in r[l]. 

In Step size, size > 1, we choose all k dimensional vectors 
with first coordinate 1 which are not linearly dependent with 
d — 2 column vectors of the constructed until now part of the 
generator matrix. Our idea is with one pass to find all proper 
vectors (all elements of U S i ze ) using all column vectors from 
the generator matrix obtained until this step using U S ize-\ (all 
proper vectors from previous step). To find all proper vectors 
Usize we take U S i ze -i and delete those which are not proper, 
with respect to the already constructed part (row ** in the 
algorithm). We find all sums of less than d — 2 columns. 
Each sum gives us one vector which is not (d — 1) -proper to 
the current step and have to be deleted. All remaining vectors 
are proper. To improve the speed of the algorithm, we pack 
each column in a computer word and use the bit operation 
XOR for computer words. The presented algorithm is much 
faster than the algorithm in ||5). 

IV. Results 

In this section, we present obtained result for codes with 
dual distance 8 and 10 using the above algorithms. With 
algorithm BRUTEFORCE we construct all codes with length 
n < 28, dimension k < 14 and dual distance at least 8, and 
all codes with dimension up to 18 and dual distance at least 
10. The summarized results for the number of inequivalent 
codes for given parameters are presented in the tables below. 
The stars in some cells mean that for the corresponding 
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parameters n and k there are codes with dual distance greater 
than the considered one. In the remaining cases, the number 
of inequivalent codes of length n and dimension k, given in 
the table, coincide with the number of optimal codes with 
minimum distance 8 (respectively 10) which have dimension 
n — k and length not larger than n. We can use the numbers 
in the tables to determine the exact number of inequivalent 
optimal codes with length n' = n and dimension k 1 = n — k 
in some of the cases. For the cells without *, the number of 
inequivalent optimal [n! = n, k' = n — k, 8] codes is equal to 
the number of inequivalent [n, k] d codes minus the number of 
inequivalent [n — 1, k — l] d codes. The calculations took about 
72 hours in contemporary PC. The number of all inequivalent 
codes with dimension 15 and dual distance 8, and dimension 
19 and dual distance 9, grows exponentially, so we could not 
calculate all cases. That is why we consider the problems for 
existence of codes with parameters [33, 18, 8] and [33, 15, 10]. 

The existence of a [33, 18, 8] code leads to the existence 
of a [33,18,8] even code and [32,17,8] even code (from 
the properties of shortened codes and Lemma |6]l and its 
dual code C32 with parameters [32, 15, c?] 8 . We know that 
£(14,8) = 28 (see Table 1) and d < 8 from the tables 
for bounds of linear codes |9 |. Theorem [8] gives us that the 
minimum distance d of C32 has to be 4 < d < 8. Using 
the algorithm EXTEND and already constructed even codes 
with parameters [28, 14, 8] 8 , [27, 14, 7] 8 , . . . , [24, 14, d > 4] 8 , 
we obtain that there are exactly two inequivalent even codes 
C\ 2 an d Cf 2 with generator matrices G\ 2 and G\ 2 and weight 



enumerators: 

1 + 124z 8 + 1152z 10 + 3584z 12 + 6016z 14 
6016z 18 + 3584z 20 + 1152z 22 + 124z 24 + z 32 . 

14 



1 + 116z 8 + 1216z 10 + 3360z 12 + 6464; 
6464z 18 + 3360z 20 + 1216z 22 + 116z 24 + z 32 . 



r 32 

V 11014z 16 
and 

h 10454z 16 



G32 — 



10000000000000000000000001111111 
01000000000000001010101101110100 
00100000000000000011100001010101 
00010000000000011011100111001110 
00001000000000010111101011111000 
00000100000000001100011101100011 
00000010000000011000111110100100 
00000001000000011000110010101011 
00000000100000011000110100010111 
00000000010000010101010100101110 
00000000001000001101100010110110 
00000000000100010001101101100011 
00000000000010000111110110110000 
00000000000001000110011010001001 

V 00000000000000100110101100101110 / 



TABLE 1 - Classification results for [n,k] d± - 8 codes 



n\k 


14 


13 


12 


11 


10 


10 














1* 


11 











1* 


4* 


12 








1* 


5* 


1 


13 





1* 


6* 


3 





14 


1* 


7* 




1 





15 


8* 


14* 


4 


1 





16 


24* 


16 


5 


1 





17 


50* 


23 


5 








18 


131 


39 


2 








19 


450 


30 


1 








20 


1863 


27 


1 








21 


11497 


13 


1 








22 


46701 


10 


1 








23 


40289 


9 


1 








24 


5177 


10 


1 








25 


536 


8 











26 


274 














27 


1 














28 


1 














29 


















G\o — 



10000000000000000000000001111111 
01000000000000001010101101100011 
00100000000000000011100001010101 
00010000000000000011010001011010 
00001000000000011000101011110100 
00000100000000011011101011111011 
00000010000000010111111110111111 
00000001000000010111110010110000 
00000000100000010111110100001100 
00000000010000011010010100101110 
00000000001000001101100010101101 
00000000000100011110101101010000 
00000000000010000110011010011101 
00000000000001001000110110010111 

V 00000000000000100110101100111001 / 



There are exactly 30481 codes with parameters [27, 18] 10 
and 11 codes with parameters [26, 18] 10 and the dual codes of 
all these codes are even. But none can be extended to a code 
with parameters [33, 19, d} 10 . This leads us to the conclusion 

Theorem 10 Codes with parameters [33, 14, 10] do not exist 
and n 2 (14, 10) = 34. 

The calculations for nonexistence of codes with parameters 
[33, 18, 8] and [33, 14, 10] took about 4 weeks in a contempo- 
rary PC. 

From Theorem 9 and Theorem 10 and tables for bounds of 
codes [9|, we have: 



None of these two codes can be extended to a code with 
parameters [33, 15] 8 . This leads to 

Theorem 9 Codes with parameters [33, 18, 8] do not exist and 
n 2 (18,8) = 34. 

It follows that codes with parameters [32, 18, 7] and 
[33, 19, 7] do not exist and the minimum distance of a putative 
[33, 19, d] 10 code has to have d < 7. From Theorem [S] and 
value of L(18, 10) = 28 we have that d has to be 5 or 6. 



Corollary 11 n 2 (19,8) = 35, n 2 (20,8) = 36, n 2 (21,8) = 
37, n 2 (22, 8) = 38, n 2 (15, 10) = 35 and n 2 (16, 10) = 36. 



References 

[1] A. Bai'g, Complexity Issues in Coding Theory. Handbook of Coding 
Theory, V. S. Pless and W. C. Huffman, Eds., Elsevier, Amsterdam, 
1998. 

[2] I. Bouyukliev, What is Q-extension?, Serdica Journal of Computing 1 
(2007) pp. 115-130. 



5 
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16 
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14 


15 
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16 











1* 
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1031* 


29 











24 


11114 


6 











25 


188572 














26 


563960 














27 


30481 














28 


11 














29 


















Algorithm EXTEND (Input: V, t): 

var U: array of sets; max,i:integer; 

found:boolean; X:set; r:array of integer; 
Procedure ext(X:set; size: integer); 
var i:integer; 
{ 

if size= t then { print(X); exit; }; 
if | U [size - 1]| = then 
{ 

if ((size>max) and (size< t)) then 
{ 

max:=size; 

if max < t then found:=true; 

}; 

exit; 

}; 

while \U[size - 1]| <> do 
{ 

if (((size + \U[size — 1]| <= max) and (max < t)) or 
((size+ | U [size — 1]|) < t) and (max — t))) then exit; 

i :— min{j : Vj £ U[size — 1]}; 

if (((size + r[i}) <= max) and (max < t)) or 
* ((size + r [i]) < t) and (max = t))) then exit; 

U[size — 1] = U[size — 
** U[size] := {vj : Vj E U[size - 1] and {XUViUVj} E P}; 

if size < t then ext({X U Vj}, size + 1); 

if found — true then exit; 

}; 



Procedure Main; 

{ 

max:=0; 

for i := |V| downto 1 do 

{ 

found:=false; 

X := {«,}; 

U[0] := {v 3 : Vj G V j > i and {X U Vj} G P }; 
ext(X, 1); 
*** r[i]:=max; 
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